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A-1 The function g(x) = f(x,0) works. Substituting 


(x,y,z) = (0,0,0) into the given functional equa- 
tion yields f(0,0) = 0, whence substituting (x,y,z) = 
(x,0,0) yields f(x,0) + f(0,x) =0. Finally, substi- 
tuting (x,y,z) = (x,y,0) yields f(x,y) = —f(y,0) — 
f(0,x) = g(x) — 80). 

Remark: A similar argument shows that the possible 
functions g are precisely those of the form f(x,0) +c 
for some c. 


A-2 Barbara wins using one of the following strategies. 


First solution: Pair each entry of the first row with the 
entry directly below it in the second row. If Alan ever 
writes a number in one of the first two rows, Barbara 
writes the same number in the other entry in the pair. If 
Alan writes a number anywhere other than the first two 
rows, Barbara does likewise. At the end, the resulting 
matrix will have two identical rows, so its determinant 
will be zero. 


Second solution: (by Manjul Bhargava) Whenever 
Alan writes a number x in an entry in some row, Bar- 
bara writes —x in some other entry in the same row. At 
the end, the resulting matrix will have all rows summing 
to zero, so it cannot have full rank. 


A-3 We first prove that the process stops. Note first that 


the product a; ---a, remains constant, because ajay = 
gced(aj,ax)lem(aj,a,). Moreover, the last number in 
the sequence can never decrease, because it is always 
replaced by its least common multiple with another 
number. Since it is bounded above (by the product of 
all of the numbers), the last number must eventually 
reach its maximum value, after which it remains con- 
stant throughout. After this happens, the next-to-last 
number will never decrease, so it eventually becomes 
constant, and so on. After finitely many steps, all of the 
numbers will achieve their final values, so no more steps 
will be possible. This only happens when a; divides a, 
for all pairs j < k. 


We next check that there is only one possible final se- 
quence. For p a prime and m a nonnegative integer, we 
claim that the number of integers in the list divisible 
by p” never changes. To see this, suppose we replace 
aj,ax by ged(aj,ax),lcm(a;,a;). If neither of aj,ax is 
divisible by p”, then neither of gcd(a;,ax),lem(a;, ax) 
is either. If exactly one a;,a, is divisible by p”, then 
Icm(a;,a,) is divisible by p” but gcd(a;,a,) is not. 
gced(aj,ax),lem(a;,a,) are as well. 

If we started out with exactly 4 numbers not divisible 
by p”, then in the final sequence a/,...,a/,, the num- 


sen? 


bers a), ,---,d;, are divisible by p’” while the numbers 


a',...,aj, are not. Repeating this argument for each pair 
(p,m) such that p™ divides the initial product a),...,dn, 
we can determine the exact prime factorization of each 
of a\,...,a),. This proves that the final sequence is 
unique. 


Remark: (by David Savitt and Noam Elkies) Here are 
two other ways to prove the termination. One is to ob- 
serve that |]; at is strictly increasing at each step, and 
bounded above by (a;---a,)". The other is to notice 
that a; is nonincreasing but always positive, so even- 
tually becomes constant; then az is nonincreasing but 
always positive, and so on. 


Reinterpretation: For each p, consider the sequence 
consisting of the exponents of p in the prime factoriza- 
tions of a),...,d,. At each step, we pick two positions 
i and j such that the exponents of some prime p are in 
the wrong order at positions 7 and 7. We then sort these 
two position into the correct order for every prime p 
simultaneously. 


It is clear that this can only terminate with all se- 
quences being sorted into the correct order. We must 
still check that the process terminates; however, since 
all but finitely many of the exponent sequences consist 
of all zeroes, and each step makes a nontrivial switch 
in at least one of the other exponent sequences, it is 
enough to check the case of a single exponent sequence. 
This can be done as in the first solution. 


Remark: Abhinav Kumar suggests the following proof 
that the process always terminates in at most (5) steps. 
(This is a variant of the worst-case analysis of the bub- 
ble sort algorithm.) 


Consider the number of pairs (k,/) with 1<k<Il<n 
such that a, does not divide a; (call these bad pairs). 
At each step, we find one bad pair (i,j) and eliminate 
it, and we do not touch any pairs that do not involve 
either i or j. Ifi< k < j, then neither of the pairs (i,k) 
and (k, j) can become bad, because a; is replaced by a 
divisor of itself, while a; is replaced by a multiple of 
itself. If k < i, then (k,i) can only become a bad pair if 
ax divided a; but not aj, in which case (k, j) stops being 
bad. Similarly, if k > j, then (i,k) and (j,k) either stay 
the same or switch status. Hence the number of bad 
pairs goes down by at least 1 each time; since it is at 
most (5) to begin with, this is an upper bound for the 
number of steps. 


Remark: This problem is closely related to the clas- 
sification theorem for finite abelian groups. Namely, 
if aj,...,d, and a\,...,ai, are the sequences obtained 
at two different steps in the process, then the abelian 


groups Z/a,;Z x--- x Z/a,Z and Z/a\Zx-+-- x Z/a,Z 
are isomorphic. The final sequence gives a canonical 
presentation of this group; the terms of this sequence 
are called the elementary divisors or invariant factors 
of the group. 


Remark: (by Tom Belulovich) A Jattice is a partially 
ordered set L in which for any two x,y € L, there is a 
unique minimal element z with z > x and z > y, called 
the join and denoted x A y, and there is a unique max- 
imal element z with z < x and z < y, called the meet 
and denoted x V y. In terms of a lattice L, one can pose 
the following generalization of the given problem. Start 
with @,...,d, €L. If i < j but a; £ aj, it is permit- 
ted to replace a;,a; by a; V aj,a; /\aj;, respectively. The 
same argument as above shows that this always termi- 
nates in at most (5) steps. The question is, under what 
conditions on the lattice L is the final sequence uniquely 
determined by the initial sequence? 


It turns out that this holds if and only if L is distributive, 
i.e., for any x,y,z € L, 
(x Az). 


A(yVz)=(xAy)V 


(This is equivalent to the same axiom with the oper- 
ations interchanged.) For example, if L is a Boolean 
algebra, i.e., the set of subsets of a given set S under in- 
clusion, then A is union, V is intersection, and the dis- 
tributive law holds. Conversely, any finite distributive 
lattice is contained in a Boolean algebra by a theorem 
of Birkhoff. The correspondence takes each x € L to 
the set of y € L such that x > y and y cannot be written 
as a join of two elements of L\ {y}. (See for instance 
Birkhoff, Lattice Theory, Amer. Math. Soc., 1967.) 


On one hand, if L is distributive, it can be shown that 
the j-th term of the final sequence is equal to the meet 
of aj, Av NGi; over all sequences 1 <i) <---<ij<n. 
For instance, this can be checked by forming the small- 
est subset L’ of L containing aj,...,@, and closed under 
meet and join, then embedding L’ into a Boolean alge- 
bra using Birkhoff’s theorem, then checking the claim 
for all Boolean algebras. It can also be checked di- 
rectly (as suggested by Nghi Nguyen) by showing that 
for j = 1,...,n, the meet of all joins of j-element sub- 
sets of a),...,@, 1S invariant at each step. 


On the other hand, a lattice fails to be distributive if 
and only if it contains five elements a,b,c,0, 1 such that 
either the only relations among them are implied by 


1>a,b,c>0 


(this lattice is sometimes called the diamond), or the 

only relations among them are implied by 
1>a>b>0, 1>c>0 

(this lattice is sometimes called the pentagon). (For a 

proof, see the Birkhoff reference given above.) For each 

of these examples, the initial sequence a, b,c fails to de- 

termine the final sequence; for the diamond, we can end 


A-4 


A-5 


up with 0, *, 1 for any of * = a,b,c, whereas for the pen- 
tagon we can end up with 0, x, 1 for any of * = a,b. 


Consequently, the final sequence is determined by the 
initial sequence if and only if L is distributive. 


The sum diverges. From the definition, f(x) =x on 
[1,e], xInx on (e,e?], xInxInInx on (e°,e], and so 
forth. It follows that on [1,c9), i is positive, continu- 


ous, and increasing. es pea Fi) y if it converges, is 
bounded below by [7° a ; it suffices to prove that the 


integral diverges. 


Write In'x = Inx and Inéx = In(In*'x) 
for k > 2; similarly write exp'x = e* and 
expkx = e's If we write y = In*x, then 


x= expty and dx = (exp*y)(exp!!y) «+ (exp! y)dy = 
x(In! x) --- (In! x)dy. Now on [exp*"!1,exp* 1], we 


have f(x) = x(In!x)---(In*!x), and thus substituting 
y=In*x yields 
exp* | 
he k-1y i s=f ae 
=52 pe 1 
k=1 k-14 io 


It follows that fy” iol on 
desired. 


=) diverges, as 


Form the polynomial P(z) = f(z) +ig(z) with complex 
coefficients. It suffices to prove that P has degree at 
least n — 1, as then one of f, g must have degree at least 
n—-1. 

By replacing P(z) with aP(z) + b for suitable a,b € 
C, we can force the regular n-gon to have vertices 
Cn, 67,---,6” for $, = exp(2zi/n). It thus suffices to 
check that there cannot exist a polynomial P(z) of de- 
gree at most n —2 such that P(i) = C! fori=1,...,n 


We will prove more generally that for any complex 
number ¢ ¢ {0,1}, and any integer m > 1, any polyno- 
mial Q(z) for which Q(i) =t! for i=1,...,mhas degree 
at least m — 1. There are several ways to do this. 


First solution: If Q(z) has degree d and leading coeffi- 
cient c, then R(z) = Q(z +1) —tQ(z) has degree d and 
leading coefficient (1—1)c. However, by hypothesis, 
R(z) has the distinct roots 1,2,...,m—1, so we must 
haved >m-—1. 


Second solution: We proceed by induction on m. For 
the base case m = 1, we have Q(1) =t! £0, so Q must 
be nonzero, and so its degree is at least 0. Given the 
assertion for m—1, if Q(i) =t! for i=1,...,m, then 
the polynomial R(z) = (t — 1)~!(Q(z+ 1) — Q(z)) has 
degree one less than that of Q, and satisfies R(i) = ft! 
fori=1,...,m—1. Since R must have degree at least 
m— 2 by the induction hypothesis, Q must have degree 
at least m— 1. 


Third solution: We use the method of finite differ- 
ences (as in the second solution) but without induction. 


Namely, the (m— 1)-st finite difference of P evaluated 
at 1 equals 


m-1 


y ay (" "acm 3) = 11-1"! £0, 


j=0 


which is impossible if Q has degree less than m— 1. 


Remark: One can also establish the claim by comput- 
ing a Vandermonde-type determinant, or by using the 
Lagrange interpolation formula to compute the leading 
coefficient of Q. 


A-6 For notational convenience, we will interpret the prob- 


lem as allowing the empty subsequence, whose product 
is the identity element of the group. To solve the prob- 
lem in the interpretation where the empty subsequence 
is not allowed, simply append the identity element to 
the sequence given by one of the following solutions. 


First solution: Put 7 = |G]. We will say that a sequence 
S produces an element g € G if g occurs as the product 
of some subsequence of S. Let H be the set of elements 
produced by the sequence S. 


Start with S equal to the empty sequence. If at any point 
the set H-'H = {hyhz : hy! ,h2 € H} fails to be all of G, 
extend S by appending an element g of G not in H~'H. 
Then Hg MH must be empty, otherwise there would be 
an equation of the form hyg = hz with hy,h2 € G, or 
g= hy ho, a contradiction. Thus we can extend S by 
one element and double the size of H. 


After k < logyn steps, we must obtain a sequence S$ = 
a,...,ax for which H~'H = G. Then the sequence 
Apogee saa sath ...,@% produces all of G and has length 
at most (2/1n2) Inn. 


Second solution: 
Put m = |H|. We will show that we can append one 
element g to S so that the resulting sequence of k+ 1 


elements will produce at least 2m —m7/n elements of 
G. To see this, we compute 


Y |HUHe| = ¥ (H+ |He|— | NHel) 
gcG gcG 


=2mn—)" HOA g| 

gcG 
=2mn—|{(g,h) € G@ he HNHg}| 
= 2mn y {geG:heHg3}| 

heH 
= 2mn — y Hh 

heH 
=2mn —m’. 


By the pigeonhole principle, we have |H UH g| > 2m— 
m~ /n for some choice of g, as claimed. 
In other words, by extending the sequence by one ele- 


ment, we can replace the ratio s = 1 —m/n (i.e., the frac- 
tion of elements of G not generated by S) by a quantity 


ee 


no greater than 
1—(2m—m? /n)/n =s?. 


We start out with k = 0 and s = 1 — 1/n; after k steps, 


we have s < (1 — 1/n)™. It is enough to prove that for 
some c > 0, we can always find an integer k < clnn 


such that 
ok 
1 1 
(1 _ ~) <-, 
n n 


as then we have n—m < 1 and hence H =G. 


To obtain this last inequality, put 
k = |2logyn| < (2/1n2)Inn, 
so that 2+! > n2. From the facts that Inn < In2+ (n— 


2)/2 <n/2 and In(1—1/n) < —1/n for all n > 2, we 
have 


1 n n 
2*in(1 = 1 
n( -) < mn 5) <—Inn, 
yielding the desired inequality. 


Remark: An alternate approach in the second solution 
is to distinguish betwen the cases of H small (i.e., m < 
n'/2, in which case m can be replaced by a value no 
less than 2m — 1) and H large. This strategy is used ina 
number of recent results of Bourgain, Tao, Helfgott, and 
others on small doubling or small tripling of subsets of 


finite groups. 


In the second solution, if we avoid the rather weak in- 
equality Inn < n/2, we instead get sequences of length 
log, (nInn) = log,(n) + log, (Inn). This is close to op- 
timal: one cannot use fewer than log, n terms because 
the number of subsequences must be at least n. 


There are at most two such points. For example, 
the points (0,0) and (1,0) lie on a circle with center 
(1/2,x) for any real number x, not necessarily rational. 


On the other hand, suppose P = (a,b),Q = (c,d),R = 
(e, f) are three rational points that lie on a circle. The 
midpoint M of the side PQ is ((a+c)/2,(b+d)/2), 
which is again rational. Moreover, the slope of the line 
PQis (d—b)/(c—a), so the slope of the line through M 
perpendicular to PQ is (a—c)/(b—d), which is rational 
or infinite. 


Similarly, if N is the midpoint of QR, then N is a rational 
point and the line through N perpendicular to QR has 
rational slope. The center of the circle lies on both of 
these lines, so its coordinates (g,h) satisfy two linear 
equations with rational coefficients, say Ag + Bh = C 
and Dg + Eh =F. Moreover, these equations have a 
unique solution. That solution must then be 


g = (CE —BD)/(AE —BD) 
h = (AF — BC) /(AE — BD) 


B-3 


(by elementary algebra, or Cramer’s rule), so the center 
of the circle is rational. This proves the desired result. 


Remark: The above solution is deliberately more ver- 
bose than is really necessary. A shorter way to say this 
is that any two distinct rational points determine a ra- 
tional line (a line of the form ax+ by+c =0 with a,b,c 
rational), while any two nonparallel rational lines inter- 
sect at a rational point. A similar statement holds with 
the rational numbers replaced by any field. 


Remark: A more explicit argument is to show that 
the equation of the circle through the rational points 


(1,91), (2,92), (%3,3) is 


tn yd 
0 =det | 3173 72 2! 
x3+y3 x3 y3 | 
Gye ees 


which has the form a(x* + y?) +dx+ey+ f =0 
for a,d,e,f rational. The center of this circle is 
(—d/(2a),—e/(2a)), which is again a rational point. 


We claim that F,(x) = (Inx —a,)x"/n!, where a, = 
Vie 1/k. Indeed, temporarily write G,(x) = (Inx — 
ay)x"/n! for x > 0 and n > 1; then lim, .9G,(x) = 0 
and Gi,(x) = (Inx—ay,+1/n)x""!/(n—1)! = Gy_i(x), 
and the claim follows by the Fundamental Theorem of 
Calculus and induction on n. 


Given the claim, we have F,,(1) = —a,/n! and so we 
need to evaluate — limy_5.. ae But since the function 
1/x is strictly decreasing for x positive, Yt, 1/k =ay,— 
1 is bounded below by f;'dx/x = Inn — In2 and above 
by fj'dx/x = Inn. It follows that lim,,.. #4 = 1, and 
the desired limit is —1. 


The largest possible radius is ee It will be convenient 


to solve the problem for a hypercube of side length 2 
instead, in which case we are trying to show that the 
largest radius is J/2. 


Choose coordinates so that the interior of the hypercube 
is the set H = [—1,1]* in R*. Let C be a circle centered 
at the point P. Then C is contained both in H and its 
reflection across P; these intersect in a rectangular par- 
alellepiped each of whose pairs of opposite faces are at 
most 2 unit apart. Consequently, if we translate C so 
that its center moves to the point O = (0,0,0,0) at the 
center of H, then it remains entirely inside H. 


This means that the answer we seek equals the largest 
possible radius of a circle C contained in H and 
centered at O. Let vj = (viq,...,Vi4) and v2 = 
(v21,---,¥24) be two points on C lying on radii perpen- 
dicular to each other. Then the points of the circle can 
be expressed as vj cos 9+ v2 sin @ for 0 < 0 < 27. Then 
C lies in H if and only if for each i, we have 
|v1,cos 8 + v2; sin O| < 1 (0<0<2n). 

In geometric terms, the vector (v1;,v2;) in R* has 
dot product at most 1 with every unit vector. Since 


4 


this holds for the unit vector in the same direction as 
(v1;, V2), We must have 
ON ig 
Vii + V9; — 1 


Gina A 


Conversely, if this holds, then the Cauchy-Schwarz in- 
equality and the above analysis imply that C lies in H. 


If r is the radius of C, then 
2 . 2 2 
2r = vit Yo v5; 
i=l i=l 


4 
= Lei +V3;) 


so r < V2. Since this is achieved by the circle through 
(1,1,0,0) and (0,0,1,1), it is the desired maximum. 


Remark: One may similarly ask for the radius of the 
largest k-dimensional ball inside an n-dimensional unit 
hypercube; the given problem is the case (n,k) = (4,2). 
Daniel Kane gives the following argument to show that 
the maximum radius in this case is se . (Thanks for 
Noam Elkies for passing this along.) 


We again scale up by a factor of 2, so that we are trying 
to show that the maximum radius r of a k-dimensional 
ball contained in the hypercube [—1,1]” is \/. Again, 
there is no loss of generality in centering the ball at the 
origin. Let T : Ré > R” be a similitude carrying the 
unit ball to this embedded k-ball. Then there exists a 
vector v; € R* such that for e} ,+++,@n the standard basis 
of R", x-v; = T(x) -e; for all x € R*. The condition of 
the problem is equivalent to requiring |v;| < 1 for all i, 
while the radius r of the embedded ball is determined 
by the fact that for all x € R‘, 


Pee): =T)2TG)= Pv 


Let M be the matrix with columns vj,...,v%; then 
MM! =r, for I, the kx k identity matrix. We then 
have 


kr? = Trace(r7I,) = Trace(MM" ) 
n 
= Trace(M’M) = y lv;|7 
i=l 
<n 


— >) 


yielding the upper bound r < “fb : 


To show that this bound is optimal, it is enough to show 
that one can find an orthogonal projection of R” onto R* 
so that the projections of the e; all have the same norm 
(one can then rescale to get the desired configuration 
of v1,.-.,V¥n). We construct such a configuration by a 
“smoothing” argument. Startw with any projection. Let 


W1,--+,Wn be the projections of €1,...,en. If the desired 
condition is not achieved, we can choose /, j such that 


[wil? < — (wil? + “+ |wal) < hwyl?. 


By precomposing with a suitable rotation that fixes 
en for h i,j, we can vary |w;|,|w,;| without varying 
|wi|? + |w;|? or |wal for h 4 i,j. We can thus choose 
such a pougnien to force one of |w;|*,|w;|? to become 
equal to +(|w |? +--+ +|wn|*). Repeating at most n — 1 
times gives the desired configuration. 


B—4 We use the identity given by Taylor’s theorem: 


h(xt+y) = » u 


In this expression, h‘?(x)/i! is a polynomial in x with 
integer coefficients, so its value at an integer x is an 
integer. 


For x =0,...,p—1, we deduce that 
h(x+p) =h(x)+ph'(x) (mod p’). 


(This can also be deduced more directly using the bino- 
mial theorem.) Since we assumed h(x) and h(x + p) 
are distinct modulo p*, we conclude that h'(x) 4 0 
(mod p). Since h’ is a polynomial with integer coef- 
ficients, we have h(x) = h'(x-+mp) (mod p) for any 
integer m, and so h'(x) #0 (mod p) for all integers x. 


Now for x =0,...,p?— 1 and y =0,..., p — 1, we write 


h(x-+ yp?) =h(x) + p’yh'(x) (mod p?). 


Thus h(x),h(x+p?),...,h(x+(p—1)p”) run over all of 
the residue classes modulo p* congruent to h(x) mod- 
ulo p?. Since the h(x) themselves cover all the residue 
classes modulo p*, this proves that h(0),...,4(p* — 1) 
are distinct modulo p?. 


Remark: More generally, the same proof shows that 
for any integers d,e > 1, h permutes the residue classes 
modulo p¢ if and only if it permutes the residue classes 
modulo p*. The argument used in the proof is related 
to a general result in number theory known as Hensel’s 
lemma. 


B-5 The functions f(x) =x-+n and f(x) =—x-+n for any 


integer 7 clearly satisfy the condition of the problem; 
we claim that these are the only possible f. 

Let g = a/b be any rational number with gcd(a,b) = 1 
and b > 0. For n any positive integer, we have 


era eas 1B) yy (S41) nbs (2) 


1 
bn 


is an integer by the property of f. Since f is differ- 
entiable at a/b, the left hand side has a limit. It fol- 
lows that for sufficiently large n, both sides must be 


equal to some integer c = f’($): f(4¢ tS) + 

Now c cannot be 0, since otherwise f(“4*~) = f(¢) for 
sufficiently large n has denominator b orhet han bn. 
Similarly, |c| cannot be greater than 1: otherwise if we 
take n = k|c| for k a sufficiently large positive integer, 
then f($) + 7 has denominator bk, contradicting the 
fact that f( antl) has denominator bn. It follows that 
e= f'(})=1. 

Thus the derivative of f at any rational number is +1. 
Since f is continuously differentiable, we conclude that 
f' (x) = 1 for all real x or f’ (x) = —1 for all real x. Since 
f (0) must be an integer (a rational number with denom- 
inator 1), f(x) =x-+n or f(x) = —x-+n for some integer 
n. 


*) = 
nl 


Remark: After showing that f’(q) is an integer for each 
q, one can instead argue that f’ is a continuous function 
from the rationals to the integers, so must be constant. 
One can then write f(x) = ax +b and check that b € Z 
by evaluation at a = 0, and that a = +1 by evaluation at 
x=I1/a. 


B-6 In all solutions, let F;, ; be the number of k-limited per- 


mutations of {1,...,n}. 


First solution: (by Jacob Tsimerman) Note that any 
permutation is k-limited if and only if its inverse is k- 
limited. Consequently, the number of k-limited per- 
mutations of {1,...,n} is the same as the number of 
k-limited involutions (permutations equal to their in- 
verses) of {1,...,n}. 


We use the following fact several times: the number of 
involutions of {1,...,n} is odd if n = 0, 1 and even oth- 
erwise. This follows from the fact that non-involutions 
come in pairs, so the number of involutions has the same 
parity as the number of permutations, namely n!. 


For n <k +1, all involutions are k-limited. By the pre- 
vious paragraph, F,,, is odd for n = 0,1 and even for 
n=2,...,k+1. 


For n > k+ 1, group the k-limited involutions into 
classes based on their actions on k+2,...,n. Note 
that for C a class and o € C, the set of elements of 
A = {1,...,4+1} which map into A under o depends 
only on C, not on o. Call this set S(C); then the size of 
C is exactly the number of involutions of S(C). Conse- 
quently, |C| is even unless S(C) has at most one element. 
However, the element | cannot map out of A because we 
are looking at k-limited involutions. Hence if S(C) has 
one element and o € C, we must have o(1) = 1. Since 
o is k-limited and o(2) cannot belong to A, we must 
have o(2) =k+2. By induction, for i = 3,...,k+1, 
we must have o(i) =k +i. 


If n < 2k+1, this shows that no class C of odd cardi- 
nality can exist, so F,, must be even. If n > 2k +1, 
the classes of odd cardinality are in bijection with k- 
limited involutions of {2k+2,...,n}, so F,% has the 
same parity as F,2,-1,x. By induction on n, we deduce 
the desired result. 


Second solution: (by Yufei Zhao) Let M,,, be the n x n 
matrix with 


1 ji-j|<k 
Mn Gi a 
(Mnt)ii {0 otherwise. 
Write det(M,,,) as the sum over permutations o of 
{1,...,n} of (Mnk)1o(1)+** Mnk)no(n) times the signa- 
ture of o. Then o contributes +1 to det(M,,) if o is 
k-limited and 0 otherwise. We conclude that 


det(M, x) = Fn (mod 2). 


For the rest of the solution, we interpret M,, , as a matrix 
over the field of two elements. We compute its determi- 
nant using linear algebra modulo 2. 


We first show that for n > 2k+1, 
Fuk = Fn—2n-1,e (mod 2), 


provided that we interpret Fo, = 1. We do this by com- 
puting det(M,, ,) using row and column operations. We 
will verbally describe these operations for general k, 
while illustrating with the example k = 3. 


To begin with, M,,, has the following form. 


SloOrrre RRR 
1 ce Oo oe od 
Cee el etl el =) 
Np RR RRR OO 
Vl RRR ROC O 

i — 


Ssl(o Oo ORR Fe 
SO OR RR Re 


* 


In this presentation, the first 2k-+ 1 rows and columns 
are shown explicitly; the remaining rows and columns 
are shown in a compressed format. The symbol @ in- 
dicates that the unseen entries are all zeroes, while the 
symbol ? indicates that they are not. The symbol * in 
the lower right corner represents the matrix F,—24—1,k- 
We will preserve the unseen structure of the matrix by 
only adding the first k + 1 rows or columns to any of the 
others. 


We first add row | to each of rows 2,...,k+ 1. 


Ssrcor OC CO 
S>crrocrcore 
Sierere OC OF 
a 
Ny S| SP SE KS OO 
Nye SF eS KS COCO Oo 

vu vr seseos 8S 


sSsrocooocooqoor 


* 


We next add column | to each of columns 2,...,k+ 1. 


Sjocdro°c’cecoco 
SJ/orrocq$octro.o°c 
Sjerere OCC & 
Oe ee =) 
ee ee 
ype Se SF KS COO oO 

vv wy Sess 8s 


Ssjooo ooor 


* 


For i = 2, for each of j =i+1,...,2k-+ 1 for which the 
(j,k +i)-entry is nonzero, add row i to row j. 


100000 0/0 
000010 0/0 
000001 0/0 
000001 1/0 
O1li1i10141)? 
0011011)? 
00010141)? 
000090? ?)x 
Repeat the previous step for i = 3,...,k-+ 1 in succes- 
sion. 
100000 0/0 
000010 0/0 
000001 0/0 
000000 1/0 
0111000)? 
0011000)? 
000100 0)? 
000090 6 0|x 


Repeat the two previous steps with the roles of the rows 
and columns reversed. That is, for i= 2,...,k-+ 1, for 
each of j=i+1,...,2k+1 for which the (j,k+/)-entry 
is nonzero, add row i to row j. 


Slocrod¢qcoo 
Sloro$¢qcrco°9o 
Slroceo°ocdcoo 


Sloqocqj#oooro 
*¥*I/S ooo ese s 


Slocoooor 
Slocvooro:e 
Slocoroeooeose 


We now have a block diagonal matrix in which the top 
left block is a (2k +1) x (2k +1) matrix with nonzero 
determinant (it results from reordering the rows of the 


identity matrix), the bottom right block is M,_—2—1,4, 
and the other two blocks are zero. We conclude that 


det(M,, x) = det(My_24—1,k) (mod 2), 


proving the desired congruence. 


To prove the desired result, we must now check that 
Fox, Fix are odd and Foy,..., Fox, are even. For n = 
0,...,k +1, the matrix M, ; consists of all ones, so its 
determinant is 1 if nm = 0,1 and 0 otherwise. (Alter- 
natively, we have F,; =n! forn =0,...,k +1, since 
every permutation of {1,...,2} is k-limited.) For n = 
k+2,...,2k, observe that rows k and k +1 of M,,, both 
consist of all ones, so det(M,,,) = 0 as desired. 


Third solution: (by Tom Belulovich) Define M,, as 
in the second solution. We prove det(M,,;.) is odd for 
n=0,1 (mod 2k + 1) and even otherwise, by directly 
determining whether or not M,, ; is invertible as a matrix 
over the field of two elements. 

Let r; denote row i of M,,. We first check that if n= 
2,...,2k (mod 2k+ 1), then M,,, is not invertible. In 
this case, we can find integers 0 < a < b <k such that 
n+a+b=0 (mod 2k+1). Put j =(n+a+b)/(2k+ 
1). We can then write the all-ones vector both as 


j-l 
y Tet1—at(2k+1)i 
i=0 

and as 


j-l 


y Tk+1—b+(2k+1)i 
i=0 


Hence M,,; is not invertible. 


We next check that ifn =0,1 (mod 2k+ 1), then M,, x 
is invertible. Suppose that aj,...,a, are scalars such 
that ayr) +---+dyr, 1s the zero vector. The m-th co- 
ordinate of this vector equals a,_,4 +++: +4@m+x, where 
we regard a; as zero if i ¢ {1,...,n}. By comparing 
consecutive coordinates, we obtain 

Am—k = Gm+k+1 (l<m<n). 


In particular, the a; repeat with period 2k+ 1. Taking 
m=1,...,k further yields that 


ky. = = 41 =0 
while taking m=n—k,...,n—1 yields 
An—2k = ++ = An—1-¢ = 9. 
For n =0 (mod 2k + 1), the latter can be rewritten as 
aj=---=a,=0 
whereas for n = 1 (mod 2k + 1), it can be rewritten as 


ay =" = ayy) =0. 
In either case, since we also have 


ai +--+ +441 =0 


from the (k + 1)-st coordinate, we deduce that all of the 
a; must be zero, and so M,,, must be invertible. 


Remark: The matrices M,; are examples of banded 
matrices, which occur frequently in numerical appli- 
cations of linear algebra. They are also examples of 
Toeplitz matrices. 


